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ON NODAL DOMAINS IN EUCLIDEAN BALLS 


BERNARD HELFFER AND MIKAEL PERSSON SUNDQVIST 

Abstract. A. Pleijel (1956) has proved that in the case of the Laplacian with 
Dirichlet condition, the equality in the Courant nodal theorem (Courant sharp 
situation) can only be true for a finite number of eigenvalues when the dimension 
is > 2. Recently Polterovich extended the result to the Neumann problem in 
two dimensions in the case when the boundary is piecewise analytic. A question 
coming from the theory of spectral minimal partitions has motivated the analysis 
of the cases when one has equality in Courant’s theorem. 

We identify the Courant sharp eigenvalues for the Dirichlet and the Neumann 
Laplacians in balls in d > 2. It is the first result of this type holding in any 
dimension. The corresponding result for the Dirichlet Laplacian in the disc in R^ 
was obtained by B. Helffer, T. Hoffmann-Ostenhof and S. Terracini. 


1. Introduction and main results 

We consider the problem of counting nodal domains of eigenfunctions of the self- 
adjoint realization of the Laplacian, —A in the unit ball in The “nodal domains” 
are the connected components of the zeroset of the eigenfunction in the ball. We 
consider the Dirichlet problem for d > 3 and the Neumann problem for d > 2 (the 
corresponding results for the Dirichlet problem for d = 2 was given in [7j). 

To be more precise, denoting by A„ the nth eigenvalue, our goal is to discuss the 
property of Courant sharpness of these operators, that is the existence of eigenval¬ 
ues A„ for which there exists an eigenfunction with exactly n nodal domains. We 
recall that Courant’s theorem says that the number of nodal domains, of an 

eigenfunction corresponding to A„ is bounded by n. Moreover, it has been proven 
that the number of Courant sharp cases must be finite, see m for the Dirichlet 
case and [19] for the Neumann case (in dimension 2 only and for piecewise analytic 
boundaries). The two first eigenvalues are always Courant sharp. We will prove the 
following. 

Theorem 1.1. The only Courant sharp eigenvalues for the Neumann Laplacian for 
the disc are Ai, A 2 and A 4 . 

Theorem 1.2. The only Courant sharp eigenvalues for the Dirichlet and Neumann 
Laplacians for the hall in d> 3, are Ai and A 2 . 

This analysis is motivated by the problem of spectral minimal A:-partitions, where 
one is interested in minimizing maxj Xi[Dj) over the family D = {Di, ■ ■ ■ ,Dk) of 
pairwise disjoint open sets in a domain fl, where Xi{Dj) denotes either the Dirichlet 
ground state energy (if we analyze the Dirichlet spectral partitions of an open set 
D) or the Dirichlet-Neumann ground state energy for the Laplacian in Dj with 
Neumann condition on dDj fl dLl and Dirichlet condition on the remaining part of 
dDj. There are now many results in the two-dimensional (2D) case. We refer to jlj 
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for a recent review. In higher dimensions much less is done, and we only know of 
the determination of all Courant sharp Dirichlet eigenvalues in the cube in three 
dimensions, m- We also know less about the properties of /s-minimal partitions in 
higher dimensions. This will not create too much problems below, because we will 
work with explicit nodal domains of eigenfunctions, which in spherical coordinates 
will be expressed as a product of an interval (in the radial direction) by a nodal 
domain of a spherical harmonics in 

In Sectionwe recall how one describes the spectrum of the Laplace operator. As 
a part of the analysis of the Neumann problem, we use and extend a recent result on 
the zeros of derivatives of the Bessel functions J^, saying that J(, and Jl^p have no 
common positive zeros if p > 0 and p > 1 are integers. This was proved by M. Ashu 
in his Bachelor thesis [2j. 

In Section we discuss Courant sharpness. As a first result, we use a symmetry 
argument to extend a result by Leydold f [T¥llT5] ) from to d > 4, saying that 
only the first two eigenvalues of the Laplace-Beltrami operator on are Courant 
sharp. We then continue towards the proofs of the Theorems ] 1.1 [ and 1.2 by reducing 
the number of cases that need special treatment by using what we call a twisting 
trick. In short, it says that if the eigenfunction is non-radial, and if the eigenfunction 
is zero on a set |a;| = p, p < 1, then one can consider the same eigenfunction, but 
where one makes a small rotation of the inner ball |a;| < p, breaking the necessary 
symmetry. We refer to Subsection |3.2| for the full details. This leaves two families of 
eigenfunctions to consider. In Subsection |3.3| we finish the proof of Theorem |1.2 


m 


the case of Dirichlet boundary condition, by using an interlacing property of zeros 
of Bessel functions. In Subsection ^ we finish the proof of Theorem o and the 
Neumann part of Theorem 1.2 We remark that the proof of Theorem o is quite 
close to the proof of the Dirichlet case for the disc [H Section 9]. 

In Section]^ we discuss the possible extension of a theorem by A. Pleijel, [T7] . 
The question is to determine if there exists a constant 7 < 1 such that, for any 
infinite sequence of eigenpairs (A^; 


, Uri 


lim sup 


p('Rn) 

n 


< 7- 


For the Dirichlet problem, this is indeed the case as proved in the paper of Berard- 
Meyer [3j, which establishes, in any dimension d>2 for bounded open sets in or 
d dimensional compact Riemannian manifolds, the existence of an explicit universal 
constant 7 (d) < 1 (extending [16]). This was also solved previously for the Neumann 
problem in dimension 2 m- 

Finally, in Section we establish new monotonicity properties of the function 

7(d). 

Remark 1.3. It would be interesting to consider the problem of minimal A;-partitions 
of the ball in three dimensions. In the case k = 3, it has been proved in [5] that the 
minimal 3-partition of the sphere is up to rotation determined by the intersection 
of with three half-planes crossing along the vertical axis with equal angle ^. It 
is natural to conjecture that the minimal 3-partition for the ball is up to rotation 
determined by the intersection of the ball with three half-planes crossing along the 
vertical axis with equal angle 
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2. Spectrum of the Laplace operator in the unit ball in 


We denote by —and —the Dirichlet and Neumann Laplace operators, 
respectively, in the unit ball in d>2. The Laplace operator —A can be written 
as 

d — ld 1/A ^ 

where r = |a;| is the radial variable and Agd-i is the Laplace-Beltrami operator, 
acting in 


Proposition 2.1 ([ 2 T|, Theorem 22.1 and Corollary 22 . 1 ]). Assume that d>2. The 
spectrum of — Agd-i consists of eigenvalues 


£{£ + d- 2 ), £en. 


The multiplicity of the eigenvalue £{£ + d —2) is given by 


^e,d ■= 




'£ + d-3\ 

. d-1 J’ 


which coincides with the dimension of the space of homogeneous, harmonic polyno¬ 
mials of degree £. 


This leads us to consider the Dirichlet and Neumann eigenvalues of the ordinary 
differential operator 

d-I d £{£ + d- 2) 
r dr ’ 

acting in L^(( 0 , 1 ), r'^^^ dr). 

The general solution to Cu = Xu is given by 

u{r) = cir V > 

where Jj/ and Y^, denote the Bessel functions of order u, and of first and second kind, 
respectively. The Bessel functions of the second kind are too singular at the origin 
to be considered as eigenfunctions. 

To state the next results, we introduce the function 

^f\r) = Ji^2e+d-2)ir), 

which is also denoted for simplicity. 


Proposition 2.2. The spectrum of —AP in the unit ball in d > 2, consists of 
eigenvalues 

= ^eN,mGN\{0}, 

where denotes the mth positive zero of the function Each eigenvalue has 
multiplicity ■ 

Proposition 2.3. The spectrum of —A^ in the unit ball in d > 2, consists of 
eigenvalues 

= ^GN,mGN\{0}, 

where denotes the mth positive (non-negative if £ = Q) zero of the function 
r i-A ^E^^\r). Each eigenvalue has multiplicity A^^d . 
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The only statement in these propositions that needs a proof is that of the mul¬ 
tiplicity of the eigenvalues. For the Dirichlet case the needed result is given in [251 
§15.28]. It says that the Bessel functions Ji, and do not have any common pos¬ 
itive zeros. This was conjectured by Bourget (1866), and follows from a deep result 
obtained by Siegel [ 22 ] in 1929. He proved that if r > 0 is an algebraic number, and 
p G Q, then Jj/(r) is not an algebraic number. 

The corresponding result for the Neumann problem was solved recently in the 
case d = 2 in Ashu’s Bachelor thesis, [2]. In this particular case the statement is 
that Jl and have no common positive zeros. Again, there is a deep result 

behind, given in [IHl page 217], which we will come back to in the proof of the first 
lemma below. 


Lemma 2.4. Assume that d > 2 and that £ & 'N. Then the positive zeros of the 
function are transcendental numbers. 


Proof. The functions Ki, (not to be mixed up with the modified Bessel functions) 
are introduced in [ 20 ] via the identity 

1 


Mr) = 


We express the derivative of E^^'^ in terms of these K functions. 


dr ^ 




V{y -|- 1 ) ^2 

L terms c 

r\i+d/2-i 


) km) 


V{£ + d/ 2 ) 


(i) 


+ rKM/2-M) 


( 2 . 1 ) 


Assume that r > 0 is an algebraic zero of r i-A -^E^f^{r). Then both A^^+d/ 2 -i(?’) and 
M+d/ 2 -M) are transcendental according to [lO] Theorem 6.3]. In particular they 
are non-zero. However, as noted in [201 page 217], also Mi+d/ 2 -i{r) is 

transcendental. But then £/r is transcendental by ( | 2 . 1 | ). Since I is an integer and r 
was assumed to be algebraic, this is a contradiction. □ 


Proposition 2.3 is a direct consequence of this lemma. 


Lemma 2.5. Assume that d > 2, £ ^ N and p G N \ {0}. Then the functions 
r i-A and r i-A have no common positive zeros. 

Before giving the proof, we recall some recursion formulas for the Bessel functions, 
valid for all p G K and positive r. 


JM) = -M'>') - Jv+M) 5 


J'M) = --Mr) + Ji^-i{r) , 

2u 

Jv+i{r) = — Jv{r) - Jv-M) ■ 
r 


Proof of Lemma 2.5. By (2.2 1-(]2.4|), we get the corresponding formulas for E^ 


( 2 . 2 ) 

(2.3) 

(2.4) 


[I] 

II 

1 

[I] 

1 

[I] 

+ 

IV 

o 


(2.5) 

[I] 

^ + d — 2„ 

= Mr) +-i-iir) , 

T 

£>1, 

(2.6) 

[I] 

2£ + d — A ^ „ 

= ^i-2{r) 

r 

, £>2. 

(2.7) 
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We divide the proof into different cases, and do the proof by contradiction, using 
recursion formulas and Lemma 12.41 

Case 1, £ = 0 and p = 1: 

If, for r > 0, ^Q(r) = = 0, then (2.5) with £ = 0 implies that Si(r) = 0, 

which contradicts Cauchy uniqueness. 

Case 2, £ = 0 and p >2: 

Assume that r > 0 is a zero of Sq and Sj,. As in Case 1, we find that Si(r) = 0, 
and so by (2.7), E! 2 (r) = —E!o(r). One application of (2.6) gives 


0 = E'Jr) = - 


p + d — 2 ^ 


Sp(r) +Sp_i(r) 


Next, we use (2.7) several times to reduce the right-hand side to an expression 
involving S 2 (r) and Si(r) only. After p — 2 applications we find a polynomial Q in 
the variable 1/r times S 2 (r) only, since Si(r) = 0. The highest degree term of the 
polynomial is 

p + d — 2 2p + d — 4 2p + d — 6 2p + d — {2p — 2) 
r r r r 

Since S 2 (r) = — E!o(r), we find that 

0 = Q{l/r)Eo{r) , 

where Q is a non-vanishing polynomial with rational coefficients. Since r is tran¬ 
scendental by Lemma [2^ Q{l/r) ^ 0. But So(r) 7 ^ 0 by Cauchy uniqueness, so we 
end up at a contradiction and conclude that Sq and Sj, have no common positive 
zero. 

Case 3, £>1 and p > 1: 


Again, assume that r > 0 is a zero of and , 


This means, using (2.5) 


and ( 2 . 6 ) respectively. 


£„ 


0 = S^(r) = -2^+1 (r) + -S^(r), 

n -^1 £+p + d — 2^ 

0 = ^^+p(r) =--- ^i+p{r) + ^i+p-i{r) 


( 2 . 8 ) 


We use (2.7) repeatedly, to reduce the second equation so that it involves only S^(r) 
and S^_|_i(r), with polynomial (in the variable 1/r) coefficients in front. The highest 
degree (in 1 /r) coefficient in front of S^_|_i(r) will, after p — 1 steps, become 

£ -\~ p d — 2 2£ 2p d — 4 2£ 2p d — 6 2£ 2p d — 2p 

r r r r ’ 

and once reduced, while calculating the determinant of the resulting system, this 
term will be multiplied with £lr (that is in front of E!/(r) in (2.8)), which will higher 
its degree (in 1/r) by one. No such term can occur elsewhere, and thus for the 
determinant of the system to be zero, r must solve a polynomial equation with 
rational coefficients, so r is algebraic. That contradicts Lemma 2.4 The other 


possibility is that 2^(r) = S^_|_i(r) = 0. But that would imply that 2^(r) = 2^(r) = 
0, which, again, contradicts the Cauchy uniqueness. □ 


3. COURANT SHARPNESS 

3.1. The result on We first analyze the case of the sphere and extends 

Leydold’s result to for d> 3. 
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Theorem 3.1. If d > 3, the only Courant sharp cases for the Laplace-Beltrami 
operator on correspond to the two first eigenvalues. 


In the proof we need the following version of Courant’s theorem with symmetry 
(see for example @1 Subsection 2.4]) which we also prove for the sake of completeness. 

Theorem 3.2. Given an eigenfunction which is symmetric or antisymmetric with 
respect to the antipodal map, the number of its nodal domains is not greater than 
two times the smallest labeling of the corresponding eigenvalue inside its symmetry 
space. 


Proof. We note that each eigenspace has a specific symmetry with respect to the 
antipodal map. An eigenfunction associated with the eigenvalue £{£ + d — 2) 
satisfies indeed 

ipei-Lo) = (-l)V^(w), Vw G . 

This is an immediate consequence of the fact that ipz is the restriction to of an 
homogeneous polynomial of degree £ oi d variables. 

With this in mind, we first assume that £ is odd, and hence let 'ipi be an eigen¬ 
function with minimal labeling u inside the antisymmetric space. Let us assume, to 
get a contradiction, that 

>2u + l. 

We note that by antisymmetry, p{'ipi) is even. Hence we would have actually 

> 2,11 + 2 . 

We now follow the standard proof of Courant’s theorem. Selecting {p{'ipi)/2 — 1) 
pairs of symmetric nodal domains, we can construct an antisymmetric function, 
which is orthogonal to the antisymmetric eigenspace corresponding to the r — 1 
first eigenvalues and has an energy not greater than the p-th eigenvalue. Using the 
mini-max characterization of the p:th eigenvalue, we get that this function is an 
antisymmetric eigenfunction which vanishes in the two remaining nodal domains. 
This gives the contradiction using the unique continuation principle. 

Next, assume that £ is even and that 'ipi is an eigenfunction with minimal labeling u 
inside the symmetric space. We assume, again to get a contradiction, that 

>2u + 1. 

We have 

pi'ipi) = p' + 2p" 

where p' is the number of nodal domains which are symmetric and p" is the number 
of pairs of nodal domains which are exchanged by symmetry. 

If p' = 0, the proof is identical to the antisymmetric one. II p' > 1, we can 
select p' — I symmetric nodal domains and p” pairs of nodal domains exchanged by 
symmetry and construct a symmetric function which is orthogonal to the symmetric 
eigenspace corresponding to the u — 1 first eigenvalues and has an energy not greater 
than the p:th eigenvalue. Here we have used our assumption by contradiction to get 
that p' — 1 + 2p" > V. We get a contradiction just as before. □ 


Proof of Theorem \3.1[ This is just an adaptation of Leydold’s proof ( [T¥l ITS] ). 

We consider the (smallest) labeling of the eigenvalue £{£ + d — 2), i.e. the smallest 
n such that A„ = ^(^ -|- d — 2). According to Proposition |2.1[ the smallest labeling 
of the eigenvalue £{£ + d — 2) is obtained by 1 if .£ = 0, 2 if .£ = 1, and 


1 + 


f£ + d-2\ 

V d-1 ) 


+ 


+ d - 3A 

V d-1 J’ 


V£>2. 
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Using Theorem 3.2 on the eigenvalue l{l + d — 2), we get 

f £ d — 3'' 


lj,{ui) < 2 


d — 1 


+ 1 


To compute this labeling we have used that for a given £, the labeling is obtained 
by adding 1 to the sum of the multiplicity associated with the I' < I with the same 
parity as 1. 

Hence, we have to check that l> 2 and d> 3, then 


d — 3 
d — 1 

_ —3 


+ 1 


< 1 + 


£ -\- d — 2 
d-l 


+ 


£ -\- d — 3 
d — 1 


(3.1) 


l+d- 
d- 

satisfied when £> 2 and d> 3 . 


Since inequality (3.1) reads 1 < ’ which is 

□ 


3.2. Twisting trick. 

Lemma 3.3. If £ > 1 and m>2 then neither nor Xf^ can he Courant sharp. 

Because the theory of minimal partitions has not been developed to the same 
extend when d > 3, we explain how the proof goes, without referring to [Jill] which 
are mainly devoted to the case when the dimension is 2 or 3. The proof below is 
somewhat reminiscent of a proof written in collaboration with T. Hoffmann-Ostenhof 
(2005), which was never published but is mentioned in [J|. 


Proof. We start with the Dirichlet situation, and omit the D in the notation. All 
eigenvalues occurring are Dirichlet eigenvalues of the Laplace operator. The domain 
will differ, and we will be explicit about that. 

Assume that we have a Courant sharp eigenvalue An = X^^m , with £ > 1 and 
m > 2. We will construct a partition of n non-intersecting open sets {Dj}^^^ in 
the ball, such that 

max Xi{Dj) < A„ . 
j 

This leads to a contradiction by the minimax characterization of the nth eigenvalue. 

Since we assume that A„ is Courant sharp, there exists an eigenfunction having 
exactly n nodal components. Moreover, this cannot be radial (since £ > 1). So 
we have ^'(r, w) = where 'ifi is a spherical harmonic. We let p\ be the 

first zero of ui^rn in ( 0 ) 1) (which exists since m > 2) and p 2 be the second zero, if 
it exists, and P 2 = 1 otherwise. The ball is naturally divided into the parts \x I < Pi 
and Pi < I a; I < 1. Next, we define the function as 


^'(a:) = 


^(Ra;), 

^(a:), 


|a;| < Pi , 

Pi < |a;| < 1 


Here R is a small rotation, constructed in such a way that the symmetry is broken. 

Let us denote hyV = Uj=i ibe twisted partition of nodal domains correspond¬ 
ing to 

We now consider a pair of nodal domains of in the form (after relabeling) 
Di := (0, pi) X n and D 2 := (pi,P 2 ) x H. The twisting leads to the pair (see 
Figure 


3.1 


middle subfigure) Di := (0, pi) x Rfl and D 2 = D 2 := (pi,P 2 ) x D. Their 
boundary is {pi} x (fl fl Rfl). The sets D\ and D 2 cannot be the 2-partition of a 
second eigenfunction in Di ^2 := Int^Di U If it was true, it would exist p G M 
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Figure 3.1. Twist showcase. The inner dashed circle correspond to 
|a;| = Pi, and the outer one to |a;| = p 2 . Left: The partition of nodal 
domains of Middle: The nodal partition of . The common 
boundary of D\ and D 2 is marked thicker. Right: The two nodal 
domains D[ and D 2 of the second eigenfunction in Di ^2 ■ 


such that jji'tJji = o R in n n Rfl. But this will imply = 'ip£ o R on by 
analyticity. We get a contradiction at the boundary of or of 

Thus, A 2 ( 1 ^ 1 , 2 ) ^ By looking at the nodal set of a second eigenfunction rxip 
in Di^ 2 , we get two new sets D[ and D 2 (the two nodal domains of 'Uip) such that 

Ai(T>i) =Ai(T>^) < A„. 

We recall that the remaining n — 2 components of the partition V have ground state 
energy A„ . This is illustrated in Figure |3T] to the right. If n = 2, then we are done. 
Below we assume that n > 2. 

We continue, by considering D[ or D' 2 , and one of its neighbors, having a boundary 
in common. Let us, for a while, denote this pair by D” and D'^. It is possible, using 
the Hadamard formula (see [12] ) to change the common boundary of and D 2 in 
such a way that two new domains D" and D 2 constructed, with 

Ai(L>'/)<Ai(5'/)<Ai(5")<Ai(L>")- 

In particular, 

max(Ai(5'/),Ai(5")) < max(Ai(L»i'), Ai(L>^')) = An • 

At this point we have constructed three domains inside the ball, with ground state 
energy strictly less than An. If n = 3, we are done. If n > 4, we continue this 
procedure recursively until all the remaining domains in the partition V have been 
modified, and find in the end a new partition V of the ball, consisting of n pairwise 
disjoint sets Dj, such that Xi{Dj) < An for all j. 

The proof in the Neumann case is unchanged. One can do the necessary defor¬ 
mations in the boundaries where the Dirichlet condition is imposed. □ 


Remark 3.4. The proof of Lemma 3^ is easier in dimension two, since we can refer 
to the “equal angle meeting” property, which is satisfied by any regular minimal 
partition as established in [7]. 

3.3. Remaining eigenvalues, Dirichlet case. 

Lemma 3.5. Let d> 3 . Then 


\D ^ \D 

^1,1 ^ "^0,2 ■ 






ON NODAL DOMAINS IN EUCLIDEAN BALLS 


9 


Proof. Denote by the mth positive zero of the Bessel function Jj,. The inter¬ 
lacing of zeros of Jj/ (see [25l §15.22]), 

0 < jv,l < ju+1,1 < jv,2 < jv+l,2 < iiy,3 < ' ' ' Vu > —1 , 

implies, with v = d/2 , that ai^i < ao ,2 j and so < 

Hence only can be (and is!) Courant sharp in the list A^^. For the sequence 
, one can use what we have proven for the sphere. Only A^^ and Af^^ can be 
Courant sharp. This completes the proof of Theorem ] 1.2 [ for the Dirichlet problem. 

3.4. Remaining eigenvalues, Neumann case. 

Lemma 3.6. Let d>2. Then 

^ 1,1 ^ ^ 0,2 • 


Proof. We show that /3i^i < / 3 o ,2 • We recall that /3i^i is the first positive zero of 
and ^ 0,2 is the first positive zero of Sq. But, according to (2.5), SQ(r) = —Si(r). 
Now, Si(r) = 2^'^/^/r(d/2 -|- l)r + 0{r^) as r —>■ 0, so, in particular Si(0) = 0 for 
all d > 2. It follows that < ^ 0,2 by the mean value theorem. □ 

As a result, cannot be Courant sharp if m > 2 . Indeed, since the eigen¬ 
functions corresponding to A^^ have precisely m nodal domains, and the labeling 
of A" is at least m -|- 1 because A^^ < A)'^^ < A ^2 • 

We continue with the eigenvalues A^. , and start with the case d = 2. We first 
see the following ordering for the eight first eigenvalues (see Figure [3^ : 

\N _ \N ^ \N _ \N _ \N ^ \N _ \N _ \N ^ ^ \N _ \N _ \N 

— "^0,1 < '^2 — ^1,1 — ^3 ^ "^4 — ^2,1 — ^5 "^6 — "^0,2 < Ay — A 3 ^ — Ag , . . . 

with corresponding number of nodal domains 

Ml = 1 ) M 2 = 2 , /X4 = 4 , /X 6 = 2 , /X7 = 6 ,... 

We observe that A ^2 < • Hence A^^ cannot have a label lower than 2.£ -|- 1 in the 

complete ordered list of eigenvalues and the corresponding eigenfunction has exactly 
2£ nodal domains. 

For d > 3 we can again use Theorem |3.1| to conclude that only the two first 
eigenvalues A^^ and can be Courant sharp. 

This finishes the proof of Theorem 1.2 in the Neumann case. 
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Figure 3.2. Nodal domains for the eigenfunctions in the Neumann 
case for d = 2. Thicker: The Courant sharp cases. The twisting 
argument excludes every case except m = 1 (row one) and £ = 0 
(column one) from being Courant sharp. 

4. On Pleijel’s Theorem 

We will discuss (the dimension-dependent) Pleijel constant 

7(d) :=C',i(Ai(Bf))-"/^<l, 

where is the Weyl constant Cd ■= {2'ir)^‘^0Jd, <jJd is the volume of the unit ball in 
and Ai(B^) is the Dirichlet ground state energy of the Laplacian in the ball B^ 
of volume 1. More explicitly, we get (see O Lemma 9]) 

_ 2‘^-^d^r{d/2f 
^ ~ ■ 

As explained in the introduction, we focus on the Neumann case. 

Theorem 4.1. For any infinite sequence of eigenpairs {Xn,Un) of the Neumann 
Laplacian in the unit ball in (d> 3), 

limsup < 7(d — 1) < 1. 

n 

We recall that 7(2) = 4/(jQ^^) and that 7(8) = ^ (see [I^). We refer the reader 
to the last section for further properties of 7(d). The Neumann case is more delicate 
but can result for the disc of the general result of Polterovich for domains in 
with piecewise analytic boundary. In m he shows that Pleijel’s theorem holds 
with the same constant as for Dirichlet, as a consequence of a fine result due to 
Toth-Zelditch [21] on the relatively small number of points at the intersection of 
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the boundary and the zeroset. To our knowledge, nothing has been established in 
dimension d> 3 for the Neumann problem. 

A natural idea is to try to control the number of nodal domains touching the 
boundary on a set with non empty interior. This was the strategy proposed by 
Pleijel [T7j for the square and more generally by I. Polterovich [19] for the 2T)-case. 
We know indeed that it is when the eigenfunction is Hence the quotient 

between the number of “boundary” nodal sets divided by the total number tends to 
0, as m —+ 00 , like 0{l/m). In this case, the “Faber-Krahn” proof works like in 
the Dirichlet case. 

Hence it remains to control the case when m < mo. In this case, £ tends to +oo 

as ^ +00. 

We know that the labeling of is larger than m x n£ where ni is the labeling 
oi £{£ + !). Hence we get 

lim sup-— < hm sup- < 7(0 — 1 ), 

where ni^rn (resp. ni) is the labeling of for the Laplacian in the ball (resp. 
of + 1) for the Laplacian on the sphere). For the last inequality, we have used 
Berard-Meyer (Pleijel like) theorem for the sphere At this stage, we have 

obtained 

lim sup < max( 7 (d), 7 (d — 1 )) < 1 . 

n 

The conclusion of the theorem is obtained using the monotonicity of 7 which will 
be established in the next section. 

5. Monotonicity of 7(d) 

We recall from [31, Lemma 9] that 7 (d) < 1. The proof of this inequality relies on 
the estimate 

r(a;) < , 

which gives the estimate for d > 18 and of numerical computations (see below) for 
d < 18. It is natural to discuss about the monotonicity of d i-A 7 (d). 

Theorem 5.1. The function d i-A 7(d) is monotonically decreasing. 

Proof. We first recall that 

Lemma 5.2. x i-A r(a;) is logarithmically convex. 


Table 5.1. Table of values of 7(d) 


d 

2 

3 

4 

5 

6 

7(d) 

0.691660 

0.455945 

0.296901 

0.192940 

0.125581 

d 

7 

8 

9 

10 

11 

7(d) 

0.081982 

0.053704 

0.035306 

0.023291 

0.015417 

d 

12 

13 

14 

15 

16 

7(d) 

0.010236 

0.006817 

0.004553 

0.003048 

0.002046 

d 

17 

18 

19 

20 

21 

7(d) 

0.001376 

0.000928 

0.000627 

0.000424 

0.000288 










12 


BERNARD HELFFER AND MIKAEL PERSSON SUNDQVIST 


We write 


and use the logarithmic convexity of the F-function, 


This implies 




d _ I 
2 ^ 


Finally, we deduce: 




{d-lf 


2{d-2) ■ 

We next estimate 3dl2-\,\l3dl2-\l2,\- A less known result by Lewis-Muldoon is 
Lemma 5.3 (|13l Formula (1.2)]). For p > 3, p i-A is convex. 

We write the convexity of ^ (for u >3) 


(5.1) 


•2 ^ ^ ( -2 _L -2 t 

— A Vjy-3/2,1 .liy-l/2,lj 


This gives 


Jv-1/2,1 


1 

< - 
- 2 


•li/-3/2,l 

Ji^-1/2,1 


+ 1 


(5.2) 


In [I], Ashbaugh and Benguria prove: 

Ud/: 


d 


d + 4 


< 


2-1,1) 


(id/ 2 , 1 )^ 


where fl is a d-dimensional domain, the right hand side inequality being attained 
for the ball and the left hand side being the Thomson inequality. In particular, for 
the d-cube, Ai(n) = d and A 2 (fl) = d -|- 3. This implies that 

3 \i/2 
d -|- 3 / 


Jd/ 2 - 1,1 


< 1 - 


Jd/2,1 


(5.3) 


To estimate the right hand side of (5.21, we use (5.3) (with d replaced by d — 1) to 
obtain 

,'2 , / „-2 


.ld/2-1,1 


1 

< - 
- 2 


.l(d-l)/2-l,l 


+ 1 < 1 - 


(5.4) 


•ld/ 2 - 1 / 2,1 3(d-i)/2,i ^ 2{d + 2) 

To our knowledge the best estimates for the zeros of Bessel functions are 

\/viv -|- 2) < jy \ < \/p + 1 (Vp + 2 -|- 1) . (5.5) 

The left estimate is available in Watson |25j, the right one was proven by Cham¬ 
bers j5j, by choosing a good trial state for the Rayleigh quotient. 

Inserting (5.1), (5.4) and the lower bound in (5.5) into the quotient 7 (d-|-l)/ 7 (d), 
we get the following bound. 


l{d + 1 ) 
l{d) 


< 2 


(d+l)2 (d-1)' 


.v/(d- l)(d + 3) d? id-2) V 2(d + 2) 


(- 


xd /2 
)/ 
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Next, we use the inequality (1 + a/xY < e“, 


(l ^ p"*' - ( 

1 1 

\ 2{d + 2)) \ 

{d + 2)l2) 


(rf+ 2)/2 


< e 


-3/4 


to get 


7(d + 1) ^ 2 


{d+lf {d-lf 


1 - 




l{d) e3/V(d-l)(d + 3) {d-2)V 2{d + 2). 

We estimate {d — l)(d + 3) > (d + 1/2)^ (valid for d > 4), and write 1 — 3/(2(d + 2)) 
as (d + l/2)/(d + 2), to find 

7(d + 1) ^ 2 (d + 1)^ (d — 1)^ (d + 2) 


7(d) ^ e3/4 d2 (d-2) (d+l/2)2 • 

We next show that the right-hand side is bounded by 2/e^/^(l -|- 5/d) if d > 4. For 
this purpose, we write 

(^+ 1 )^ (^- 1 )^ {d + 2) _ 

(P (d-2) (d +1/2)2 i 
d5(-4 - 12/d + 39/d2 + 41/d=* - 2/d4 - S/d^) 


4(d-2)d2(d+ 1/2)2 


Since 


-4-12/d + 39/d2+41/d^-2/d^-8/d® < -4 + 39/d2 + 41/d^ Vd > 1, 
and —4 -|- 39/d^ -|- 41/d^ is monotonically decreasing and equals —59/64 for d = 4, 


we find that 


Thus 


(d+l)2(d-l)2 (d+2) 


d2 (d-2) (d + 1/2)^ 


< (1 + 5/d), Vd>4. 


<^(l + 5/d), Vd>4. 


7(d) e-^ 

The numerical approximation 2/e^/^ « 0.945 implies that 2/e^/^ < 95/100, and so 


+ Vd>4 


7(d) 


100 V d 


The right-hand side is monotonically decreasing in d, and equals one for d = 95. 
Hence, 

7(d+ 1) 


7(d) 


< 1, Vd>95. 


The remaining cases, 2 < d < 94, are covered numerically using Mathematica, and 
the quotient is plotted in Figure |5.1[ This finishes the proof of Theorem ^ □ 


Remark 5.4. As shown in Figure EH a classical asymptotics for jV,! gives, after 
observing that 

lim 

d^+oo 7(dj y d-).+oo 7(dj 
the following limiting value for the quotient 7 (d + l)/ 7 (d): 


lim 

d —^~|“00 


lid + Y 

7(d) 


2 

e 
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1 -- 


2/e - 


2 20 40 60 80 94 

Figure 5.1. The quotient 'y{d + l)/'y{d), for 2 < d < 94 . The limit is 2/e . 
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